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Superconductivity on the surface of topological insulators is known to be anisotropic and 
unconventional in that the symmetry is the mixture of s-wave and nodeless p-wave component. 
In contrast to Anderson's theorem for the insensitivity of the s-wave superconducting critical 
temperature to the nonmagnetic (time-reversal symmetric (TRS)) impurities, anisotropic su- 
perconductors including nodeless p-wave one are in general fragile even with small concentration 
of the TRS impurities. By employing the Abrikosov-Gor'kov theory, we clarify that this type 
of unconventional superconductivity emergent on the surface state of the strong topological 
insulators robustly survive against TRS impurities. 

KEYWORDS: topological insulator, helical Dirac electron, unconventional superconductivity, impurity scat- 
tering, time reversal symmetry 



Topological insulator(TI) is a new quantum state of 
matter. 1 - 4 ) TIs are fully gapped in bulk as ordinary in- 
sulators but also have topologically protected conducting 
states on their boundaries. For example, two-dimensional 
(2D) TIs have one-dimensional ballistic conducting states 
and this conducting states are found in HgTc quantum 
wells. 5 ) Moreover, a class of three-dimensional (3D) TI 
namely, strong TI, 2 ) has been predicted to have metal- 
lic 2D surface states, which have recently been observed 
by angle-resolved photoemission spectroscopy 6,7 ) Such 
surface states consist of so-called helical Dirac electrons, 
which occupy a single spin state (pointing specific di- 
rection s ) depending on the electrons' momentum p as 
sex. p. 

One of the most notable properties of strong TI is the 
robustness of the metallic surface states to TRS impuri- 
ties. This robustness originates from the absence of back- 
scattering processes when the system has the time rever- 
sal symmetry. In the non-interacting electron systems, 
back-scattering processes due to impurities play impor- 
tant roles in forming localized states that do not con- 
tribute to the electric current at zero temperature. 8 ) In 
2D electron systems, such localization leads to insulating 
states at zero temperature in the thermodynamic limit. 9 ) 
However, when back-scattering processes are prohibited 
from some symmetric reasons, metallic states persist. 

In other words, weak anti-localization is observed in 
TIs when we focus on weakly localized electron systems. 
Impurity scatterings usually induce negative quantum 
corrections to conductivities. However, in systems with a 
spin-orbit coupling (SOC) such as strong TIs, the quan- 
tum correction to the conductivity becomes positive and 
results in weak anti- localization effects. 10 ) Furthermore 
when Dirac electrons, which are not necessarily helical, 
are scattered by impurities, accumulated ir Berry phase 
around Dirac cone yields weak anti-localization. 11 ) In the 
most famous Dirac electron system, graphene with four 
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Dirac cones, there are inter- valley scatterings, which lead 
to localization in the presence of disorders. In contrast to 
the graphene, helical Dirac electrons on surfaces of strong 
TIs arc composed of odd numbers of Dirac cones and they 
show robust anti-localization. 2 ) 

Recently possible anisotropic unconventional super- 
conductivity (SC) with the mixture of s-wave and node- 
less p-wave component in TIs has been studied. It is 
known that even an s-wave attractive interaction nec- 
essarily induces the pairing potential composed of the 
mixture of s-wave and p-wave component on the surface 
of TI. 12 ) Moreover by proximity effect with s-wave super- 
conductors, SC with the same symmetry above is induced 
on the surface of TI in contact with the conventional su- 
perconductors when the Fermi energy is away from Dirac 
point. 13, 14 ) 

Such unconventional SC attracts much attention from 
viewpoints of applications. The unconventional SC on the 
surface of TI is proposed to apply for quantum computa- 
tions using Majorana fermions caused by the proximity 
effect between a superconductor and the surface states of 
TI. 13 ) Introducing SC on the surfaces of TI has been a 
challenge for experimental researches. 15 ) 

A fundamental property of SC is their impurity ef- 
fects. It is well known that the s-wave SC is robust to 
the TRS impurities 16 ) while anisotropic unconventional 
SCs such as ri-wave or p-wave SC are fragile to TRS im- 
purities in general. 17 " 20 ) A simple extension also leads to 
the fragility of the nodeless p-wave as well. 

Although the unconventional SC on the surface of 
TI have s-wave pairing potential component, they nec- 
essarily have p-wave component with the same ampli- 
tude. Therefore the response to TRS impurity concentra- 
tion can be different from ordinary s-wave SC. Moreover 
due to the anti-localization which helical Dirac electrons 
show, the impurity effect on this type of SC is nontrivial. 
In this letter, in contrast to the widely accepted fragilities 
of the unconventional superconductors, we reveal that an 
unconventional SC stabilized on the surface of TI is ro- 
bust to the TRS impurity scattering. This conclusion is 
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obtained by analyzing the dependence of the critical tem- 
perature T c on the TRS impurity concentration by the 
Abrikosov-Gor'kov (AG) theory 22 ). 

In order to study how TRS impurities affect the SCs 
on the surface state of strong TI, we introduce an ef- 
fective model of the surface state. Our model is a 2D 
helical Dirac electron system with an s-wave attractive 
interaction and in the presence of small concentration of 
impurities that cause on-site TRS scatterings. Therefore 
our Hamiltonian consists of three parts i.e. a 2D helical 
Dirac electron dispersion, Ho, an s-wave attractive inter- 
action term, H- lnt and an on-site TRS impurity scattering 
term, H lmp : 

H = H + H int + H[ mp . (1) 
The 2D helical Dirac electron dispersion H is written 

as 

H = Y, ci (k)(v F (T ■ k- fil)c(k) (2) 
fc 

with the Fermi velocity vp > 0. In this paper, we set the 
Fermi energy /x > for simplicity. At /i = case, we 
have to consider two branches of Dirac electrons because 
the both branches contribute to SC. Moreover, when fj, 
is closer to 0, the SC is harder to be observed because 
the density of states at the Fermi energy gets fewer. Pauli 
matrix er describes electron's spin and c(k) — (cfe-j- Ck±) T . 
We consider the 2D Dirac Hamiltonian of the conduc- 
tion electrons with only one Dirac cone for simplicity. 
We note that a model with single Dirac cone is suffi- 
cient to understand the essential properties of the surface 
states of strong TIs. By using a unitary transformation 

4.r = (4 T + Tei6kc ki)l"^ ( t = + 01 -) > this term is 

diagonalized as 

H = Y / dHk)(v F \k\T z ~ f i)d(k), (3) 

fc 

where t z is Pauli z matrix describing branches of Dirac 
electrons and d(k) = (d k + dk-) T ■ The index r (r = ±) 
represents branches of Dirac electrons. Here the "+" ( "-" ) 
branch represents the branch above (below) Dirac point. 
In the unitary transformation, an angle parameter 9 k = 
&rg(k x +ik y ) is introduced. The operators d and d) satisfy 
the anti-commutation relation {d kT ,d k , T ,} = S kk '8 TT i. 

Then we introduce an s-wave attractive interaction 
-ffint with an energy cutoff. For example, the attractive 
interaction may originate from an electron-phonon cou- 
pling. We assume that H lri t is written as 

H int = ^ J! Vkk>ss>J_ ks ci s ,c k , s ,c- k , s . (4) 

fc.fc' ,s,s' 

In this equation, we assume that 

t/ _ T/ _ f-9(£k,£k< € [-w c ,w c ]) , , 

T/ _ T/ _ f-9' G [-Wc.Wc]) / f>\ 

yfefe'tt — V kk'U — 1 n It t h r n J W 

with the cutoff uj c <C fJ., and g,g' > 0. Here we define 
the energy as the energy of the "+" branch measured 



from the Fermi energy as = vp\k\ — /i. We note that, 
from the condition uj c <C /U, the interaction only affects 
electrons on the "+" branch. Therefore, we neglect "-" 
branch and we write d^P as d^) below. Then the inter- 
action term is written as 

H mt ~ -£j2 ei{<>h '~ eh)d -k d Uk'd- k <. (7) 

* 

The summation represents that taken in the region 

k.k' 

€k, Cfe' G [-w C) cjJ. The terms including Vfcfc'tf and Vfefe'4.4. 
vanish because they must be odd functions of k. On the 
other hand, the terms including Vkk"\--l an d ^fefe'tJ. survive 
because they are even functions of k. 

Here we introduce an on-site TRS impurity scattering 
term as 

JVi 

ffb-p = I E E e- i9 - R, J(k + Q)c(k), (8) 

i— 1 k,q 

where N\ is the number of the impurities in the system, 
S is the size of the system, and Ri (i = 1, • • • , iVj) is the 
impurity location. Because on-site impurity scatterings 
are assumed here, scattering amplitudes do not depend 
on momentum transfer q. By using the unitary transfor- 
mation introduced above eq.(3), eq.(8) leads to 

Ni 

Himp = I EE e- iq Ri P(0k - 0k +q )d{ +q dk, (9) 

i— 1 k,q 

where P{6) = e l9 / 2 cos(0/2) is a phase factor specific to 
Dirac electron systems. The phase factor P(9k — 0k+ q ) 
contributes to the 7r Berry phase and leads to anti- 
localization effect of single Dirac cone systems. n ) 

We introduce a SC paring potential and write down 
the BCS mean-field Hamiltonian. By introducing a pair 
potential 

\{k) = e- ie «gY,z l<>h '{dk>d- kl ), (10) 
fc' 

our interaction term is approximated as 

Hint ^ - E [ A ( k ) d -k4 + A*(fc)dfcd-fc] , (II) 
fc 

/ 

where the summation represents that taken in the 
fc 

momentum space such that £ k G [— w c ,w c ] and k x > 0. 
We linearize the gap equation in terms of A because 
we only consider T c . The pair potential represented by 
eq.(10) is valid when we can neglect the mixture of "+" 
and "-" branches. This pair potential depends on the mo- 
mentum k and its dependence is represented as 

A(k) = Ae- i9h . (12) 

Therefore the equation A(fe) = — A(— fc) holds. This 
pair potential resembles one appearing in the spinless chi- 
ral p-wave superconductor. 21 ) 

This SC order is unconventional because it is composed 
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of a mixture of two different pairing symmetries. The 
pairing potential matrix in terms of the original electrons 
operator c is defined from the equation 



EE 

k Si, 52 



[d k8i A slS2 (k)i 



,t 



/i.e. 



(13) 



where 

A(fe) = A s (fe) + A 4 (fe) 
A/0 1 



A s (fe) 



2 V— 10 



A t (fe) 



A 
~2 



-e lf,fc 



.(14) 



Therefore the SC induced by iJ; n t is the mixture of singlet 
s-wave pairing and triplet nodeless p-wave pairing with 
the same amplitude. 12 ) The s-wave pairing term A s has 
no k dependence while the p-wave term A t has a phase 
determined from the direction of k. 

In order to analyze T c , we construct Gor'kov equations 
and linearize the gap equation to determine T c . 22 ' Be- 
cause the long-range SC order does not develop in 2D 
systems at finite temperature, T c calculated from the 
mean-field theory provides a criteria of the Berczinskii- 
Kosterlitz-Thoulcss (BKT) transition for development of 
the quasi-long-range order. 23 ' 24 ' 

First, we introduce two thermodynamic Green's func- 
tions G fc , fc /(r) = -(T T d k (T)d{,) and F* k ,(r) = 

(T T d_ fc (r)dL), where G (F*) is the normal (anoma- 
lous) Green's function, which describes the dynamics 
of Cooper pairs. From the equations of motion for two 
Green's functions, we can construct the Gor'kov equa- 
tions as 



h.k' = {iu n - £,k)Gk,k'{iu n ) + A(k)F kk ,(iuj n ) 

G q , kl (icj n )e i ^- k > B -P(e q -e k ), (15) 



Ni 
z—l q 



and 



= (iw„ + £ k )F, 

N 



it 



EE f . 



q,k' 



k:k ,(ioJ n ) + A*(k)G k ,k'(iuJn) 
{ioj n )e^- k > Ri P(e k -9 q ), (16) 



i=l 



where u n = (2n + l)nT is the fcrmionic Matsubara fre- 
quency, and T is the temperature. Moreover, we assume 
that A(fc) has a value independent of the frequency when 
G [— cjci w c]- By linearizing in terms of A, eq.(15) leads 

to 



S k , k ' 



{iu n - ik)G k M'{iu n ) 
? E E G «,* (^n)e^- k >^P(6 q 6 k ). (17) 



i=l 



By using a perturbation series expansion with respect to 
u/S, the Green's function G is represented as 



Gfc, fc ' (iui n ) = ^2 G k n ' k , (iu) n ) 



(18) 



where 



and 



AT; 



q i—1 



xP(6 q -e k )G ( q n - 1) (icj n ) (n=l,2,- 



(19) 



The non-perturbative Green's function G^(ia; n ) = (i,Lo n — 
^k) 1 is introduced above. By substituting eq.(18) and 
eq.(19) into eq.(16) and using eq.(17), the anomalous 
Green's function is obtained as 

F k,k' ( iuj n) = ^ A* (q)G- q - k (-iu) n )G q , k < (iw n ).(20) 



On the other hand, from eq.(10), we obtain 



A * = 9 T EE e ^'^*^ 

iuj n k' 



ILJ 7 - 



(21) 



where A*(fc) = A*e l9h . From cqs.(20) and (21), we ob- 
tain the equation to determine T c as follows: 



(gTc)- 1 =^2 F (n,m) 



F(r 



EEE 1 



x G 



(n) 
-k,-k 



,{-iu n )G k m k \{iu n ) 



(22) 



Then we analyze dependence of T c on the impurity 
concentration n, = Ni/S by using the perturbative AG 
theory, where rij is the small parameter. 22 ) According 
to the AG theory, we neglect spatial correlations be- 
tween different impurities. This approximation i.e., the 
impurity average, is justified when the impurity concen- 
tration is small. By the impurity average operation, the 

1 Nl 

term — e l i-{ R i R i) j s replaced by rii, and only 0(rii) 

terms are retained. We perform the impurity average of 
the right hand side of eq.(22) for each set of (n, m), step 
by step, as follows: We only retain the terms in eq.(22) 
that satisfy n, m < 2 and n + m < 2 because these terms 
are of the lowest order in the n^-dependence of T c . The 
non-perturbative term i.e., F(0, 0), which is illustrated in 
Fig. 1(a), is calculated as, 



F(0,0) = )]^G»K)G«_ fc (- !W „). 



(23) 



The phase factor vanishes because only the term satisfy- 
ing k = k' remains. The term F(0, 1) and -F(1,0) shown 
in Fig. 1(b) are negligible for the estimation of T c be- 
cause these terms just cause a constant self-energy shift 
and do not contribute to relaxation processes due to im- 
purity scatterings. Then F(l,l), which corresponds to 



n=0 
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the diagram illustrated in Fig. 1(c), is represented as 



concentration, 



riiU 
~S~ 



££G° fc (z Wn )G° fc ,( lWll ) 



iuj n k.k' 



(24) 



The diagram illustrated in Fig. 1(d) represents F (0, 2) or 
F(2, 0), and is given by 



F(0,2) = F(2,0) 



riiU 
~S~ 



EEE 



cos 



iuj n k q 

x Gi{iuj n )Gl(iLo n )G\(-iuj n ). (25) 
To combine these terms, F(0, 0), F(l, 1), F(2, 0), and 






(d) 



Fig. 1. Diagrams for terms of the lowest order in eq.(22). (a), (b), 
(c) and (d) correspond to (n,m) = (0,0), (n,m) = (0, 1) or(l, 0), 
(n,m) = (1,1) and (n, m) = (0,2) or (2,0), respectively. Wavy 
lines represent impurity scatterings. 



F(0, 2), we obtain 



1 = ^2: 



l + (2r 2 | W „|)- 1 



iuj n k 



[l + (2r 1 K|)-i]^+e| 



(26) 



where n and T2 are relaxation times which are related 
to F(2,0) and F(l, 1), respectively. Both quantities are 
calculated as 

_t_ = Trr^^^^^^^ = ^ n . u 2 N ^ (2?) 



and 



1 irniU 2 



2t 2 



s V 



<5(C?) = ^M 2 iVo, (28) 



where iV = J2k $(£,k)/S is the density of states at Fermi 
energy. According to the AG theory, when the concen- 
tration n% is small, T c is obtained from the equation. 22 ) 

T c ( ni ) = T c (0) - — (29) 

where 'F is the digamma function. In these equations, 

2e UJc 2_ 

T c without impurities is T c (0) = e n N o and t s is 

7T 

obtained from (r,)- 1 = (2-ri)- 1 - (2T2)" 1 . Then, from 
eqs.(27) and (28), we obtain (t,) -1 =0, up to the lowest 
order of Hj , and we conclude at least for a small impurity 



T c {0) - T c {i 



(30) 



Such cancellation of the relaxation times does not occur 
in the case of <i-wave 17 ' or p-wave 19 ' 20 ) SC orders. The 
present SC order is anisotropic because it is composed of 
s-wave and p-wave order. Therefore, our result supports 
that anisotropic but robust SC orders can exist on the 
surface of TIs. 

We have shown that unconventional SCs induced by 
the s-wave attractive interaction on surfaces of TIs are 
robust to TRS disorders. This conclusion is achieved by 
calculating dependence of T c on TRS impurity concen- 
tration, where T c provides a criteria of the BKT tran- 
sition for the quasi- long-range order in 2D. Unconven- 
tional SC studied in the literature, such as the d-wave 
and chiral p-wave SC, is sensitively suppressed through 
scatterings by a tiny concentration of impurities because 
of the phase factor of the pairing potential. In marked 
contrast, the unconventional SC on the surface of TI is 
robust because of the cancellation of two phase factors, 
one from the pairing potential and the other arising when 
a Dirac electron is scattered by a TRS impurity. The ro- 
bustness of the unconventional nodeless p-wave SC may 
favorably be tested in experiments. This result is applica- 
ble to SC on TI by proximity effect 13, 14 ) as well because 
it has the same Cooper pairing symmetry as the SC con- 
sidered in this letter. An issue left for future is the case 
of anisotropic attractions, where SC order with different 
symmetries may occur. 
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